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SECTION  I 


INTRODUCTION 


1.  BACKGROUND 

In  19o8  the  United  States  Air  Force  initiated  a  comprehensive  research 
program  to  develop  a  rational  pavement  evaluation  procedure.  Three  capabili¬ 
ties  expected  from  this  evaluation  procedure  were  as  follows  (ref.  1): 

(1)  A  technique  which  could  he  used  to  determine  the  stresses, 
strains,  and  deflections  at  any  point  within  a  pavement 
structure  loaded  by  the  multiple  wheels  of  various  gear  con¬ 
figurations  of  modern  aircraft 

(2)  An  analytical  capability  by  which  the  influence  of  varying 
the  strength  and  thickness  of  the  various  pavement  compo¬ 
nents  could  be  compared  in  terms  of  the  stress  and  strain  re¬ 
sponse  of  the  pavement  system 

(3)  Applicability  of  this  procedure  to  both  rigid  and  flexible 
pavements,  composite  pavements,  and  pavement  overlays 

The  Air  Force  Weapons  Laboratory  (AFWL)  in  1969  assigned  the  task  of  de¬ 
veloping  an  analytical  model  to  the  Naval  Civil  engineering  Laboratory  (NCEL) . 
Crawford  at  NCFL  developed  a  finite  clement  computer  code,  called  Airfield 
Paoernent  (AFP A  V)  code,  which  models  the  layered  pavement  system  as  a  prismatic 
solid  (ref.  2).  The  theoretical  basis  for  this  analytical  model  is  derived 
from  Herrmann's  three-dimensional  elasticity  solution  for  continuous  beams 
(ref.  3j.  The  central  part  of  this  code  is  based  on  a  finite  element  program 
which  Herrmann  developed  for  three-dimensional  elasticity  analysis  of  period¬ 
ically  loaded  prismatic  solids. 

While  the  AFPAV  code  was  being  developed  by  NCFL,  the  Fric  H.  Wang  Civil 
engineering  Research  Facility  (CERF)  operated  by  the  University  of  New  Mexico 
for  the  U.S.  Air  Force  evaluated  some  existing  computer  codes  to  determine 
their  apnlicability  for  theoretical  analysis  of  general  pavement  systems.  The 
results  of  this  evaluation  are  reported  in  reference  4.  The  central  part  of 
the  AFPAV  code  was  compared  with  computer  codes  such  as  BISTRO  (ref.  5),  VISAB3 
(ref.  6 J ,  and  WIL67  (ref.  7).  The  AFPAV  code  was  found  to  be  the  most  suitable 
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and  efficient  of  these  codes  for  obtaining  a  complete  picture  of  the  stress 
distribution  in  layered  pavement  systems,  both  rigid  an:'  flexible,  and  for 
plotting  pavement  deflection  basins  under  multiple-wheel  aircraft  gears  (c.g., 
C-5\  and  Boeing  747).  furthermore,  the  AFPAV  code  can  be  used  to  analyze  ef¬ 
ficiently  certain  complex  pavement  structures,  which  could  be  idealized  as 
prismatic  solids,  without  having  to  use  three-dimensional  finite  element  pro¬ 
grams.  Since  the  AFPAV  code  uses  only  a  two-dimensional  finite  element  mesh, 
hut  solves  special  types  of  three-dimensional  problems,  it  is  also  referred  to 
as  an  extended  two-dimensional  code.  In  view  of  the  great  promise  of  the  AFPAV 
code  it  was  decided  to  examine  its  full  potentialities  by  using  it  for  further 
analyses  of  structural  response  of  pavement  systems  subjected  to  heavy  aircraft 
loads. 

A  linear  clastic  model  is  used  for  the  entire  pavement  system  in  the  origi¬ 
nal  version  of  the  AFPAV  code  (refs.  2,4).  Since  the  inadequacy  of  this  sim¬ 
plified  constitutive  relationship  to  represent  all  of  the  pavement  components 
is  realized,  the  code  is  currently  being  modified  to  consider  not  only  nonlinear 
material  properties  but  also  vertical  joints  in  rigid  pavement  systems  and  dis¬ 
continuities  at  layer  interfaces,  etc.  This  task  has  been  given  to  NCfiL  by  AFWL, 
and  is  expected  to  be  completed  in  FY73 .  After  checking  out  the  performance 
of  the  modified  AFPAV  code  and  establishing  its  usefulness  as  an  analytical  tool 
for  airfield  pavement  evaluation,  a  report  containing  complete  documentation  of 
the  final  version  of  the  code  along  with  a  user's  manual  and  computer  graphic 
displays  of  the  analytical  results  of  some  sample  problems  will  be  published. 

2.  OBJECTIVES 

The  first  objective  of  this  research  effort  was  to  conduct  a  theoretical 
analysis  of  airfield  pavement  systems  to  determine  the  pavement  response*  to 
multiple  wheel  loads  using  the  finite  clement  analysis  technique.  The  rigid 
and  flexible  pavements  of  the  Multi-Wheel  lleavy-Gcar  Load  (MW1IGL)  test  sections 
constructed  and  tested  at  the  Waterways  Experiment  Station  (WHS)  were  analyzed 
using  the  original  AFPAV  code.  The  predicted  responses  were  compared  with  the 
measured  stresses  and  deflections  in  the  test  sections. 

The  second  objective  was  to  perform  a  parametric  study  of  the  properties  of 

the  different  pavement  components  of  the  MW11GL  test  sections  using  the  first 

version  of  the  AFPAV  code  to  determine  the  changes  in  pavement  response  due  to 

changes  in  the  properties  of  the  different  pavement  components. 

& 

Elastic  deflections,  strains,  and  stresses. 
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While  implementing  the  above  objectives,  the  capabilities  and  the  limi¬ 
tations  of  the  original  AFPAV  code  were  examined  in  some  detail.  Its  appli¬ 
cations  to  theoretical  analysis  of  pavement  structures  were  then  assessed  by 
flj  analyzing  the  MWliGL  test  sections  and  comparing  the  computed  results  with 
the  experimental  data,  and  (2)  conducting  parametric  studies  of  both  flexible 
and  rigid  pavement  systems. 

3.  APPROACH 

The  objectives  of  this  research  effort  which  relate  only  to  finite  ele¬ 
ment  analysis  of  pavement  structures,  were  accomplished  in  two  phases  using 
the  original  version  of  the  AFPAV  code  in  which  a  linear  elastic  model  was 
used  for  the  entire  pavement  system. 

a.  Phase  I- -Evaluation  of  AFPAV  Code 

The  AFPAV  code  consists  of  three  parts:  (1)  AFPRE  (preprocessor), 

(2)  AFPAV  (main  program),  and  (3)  AFPOST  (postprocessor).  The  different  parts 
of  the  code  were  evaluated  by  analyzing  the  MWHGL  test  sections  at  WES.  Al¬ 
though  the  finite  element  analysis  technique  can  handle  nonlinear  material  prop 
erties,  only  linear  elasticity  theory  was  assumed  in  the  first  version  of  the 
AFPAV  code  which  was  used  for  the  analyses  conducted  in  this  effort.  The  modul 
of  elasticity  of  the  various  pavement  layers  were  estimated  from  the  available 
laboratory  and  field  data.  The  computed  response  (i.e.,  deflections  and  stress 
es)  was  compared  with  the  pavement  response  measured  by  WES  in  the  full-scale 
testing  of  pavement  sections  loaded  with  a  12-wheel  assembly  of  the  main  land¬ 
ing  gear  of  the  C-5A  aircraft. 

b.  Phase  II- -Parametric  Study 

The  AFPAV  code  was  used  to  perform  a  parametric  study  to  determine 
the  effects  of  varying  the  properties  of  pavement  structure  layers  on  the  size 
and  shape  of  the  deflection  basin  of  a  pavement  under  multi-wheeled  aircraft. 
Rigid  and  flexible  pavement  systems  of  the  MWHGL  test  sections  were  studied. 

The  thicknesses  of  the  pavement  layers,  as  well  as  the  elastic  constants  (i.e., 
the  modulus  of  elasticity,  E,  and  Poisson's  ratio,  v) ,  were  varied  and  the  ef¬ 
fects  of  changing  these  properties  were  examined.  In  addition  to  pavement  de¬ 
flection  basins,  variations  of  vertical  and  horizontal  normal  stresses,  and 
shear  stresses  caused  by  variations  in  pavement  layer  properties  were  examined. 


3 


SECTION  If 


DESCRIPTION  OE  AFPAV  CODE 

The  finite  element  structural  analysis  technique  used  in  the  AFPAV  code  is 
discussed  in  detail  in  reference  2  and  also  described  briefly  in  reference  4. 
Therefore,  only  a  brief  account  of  the  salient  features  of  the  theory  is  given 
here.  This  section  also  describes  the  three  component  parts  of  the  AFPAV  code 
as  well  as  the  special  features  of  the  code  to  show  how  it  can  perform  analyses 
of  certain  pavement  problems  that  cannot  be  accomplished  by  other  codes  such  as 
BISTRO  and  WIL67. 

1.  THEORY  OF  PRISMATIC  SOLID  ANALYSIS 

The  AFPAV  code  is  based  on  a  finite  element  analysis  computer  program  orig¬ 
inally  written  by  Herrmann  for  performing  three-dimensional  elasticity  analysis 
of  periodically  loaded  prismatic  solids  (ref.  8).  Before  writing  this  program, 
Herrmann  presented  the  theory  of  periodically  loaded  structures  in  a  paper  on 
three-dimensional  analysis  of  continuous  beams  (ref.  3),  in  which  he  defines 
a  prismatic  solid  (fig.  1)  as  a  body  (1)  that  is  infinite  in  extent  in  the 


Figure  1.  Prismatic  Solid  (Continuous  Beam) 


longitudinal  direction,  z;  (2)  whose  cross  section  is  identical  for  all  values 

of  z;  and  (3)  whose  material  properties  do  not  vary  in  the  z-direction.  The 

cross-sectional  dimensions,  compared  to  the  length  of  the  prismatic  solid,  21, 
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are  such  that  it  is  necessary  to  analyze  the  solid  by  the  three-dimensional 
elasticity  theory.  Furthermore,  the  analysis  is  restricted  to  those  problems 
in  which  the  spatial  dependence  of  the  loading  may  be  approximated  as  being 
periodic  in  the  z-dircction.  The  period  of  the  loading,  however,  can  be  made 
sufficiently  large  so  that  the  effects  of  isolated  single  loads  (or  groups  of 
closely  spaced  loads)  can  be  analyzed  with  reasonable  accuracy.  For  example, 
the  prismatic  solid  analysis  may  be  used  to  determine  the  effect  of  a  point 
load  on  a  long  tunnel  or  culvert  by  considering  a  series  of  point  loads  sepa¬ 
rated  by  sufficiently  long  distances,  2L,  to  prevent  or  minimize  any  interac¬ 
tion  (figs.  2  and  3). 

In  the  AFPAV  code  the  layered  pavement  system  is  idealized  as  a  prismatic 
solid  which  is  described  in  a  rectangular  Cartesian  coordinate  system  (x,  y, 
and  z);  the  z-axis  is  assumed  to  be  perpendicular  to  the  transverse  section  of 
the  pavement  structure  (fig.  4). 

The  analysis  of  the  effects  of  any  arbitrary  gear  loads  acting  on  a  pave¬ 
ment  structure  is  made  tractable  by  expressing  the  concentrated  loads,  which 
arc  functions  of  z,  as  a  Fourier  series  in  z.  Figure  5  shows  the  configuration 
of  a  12-wheel  assembly  of  a  C-5A  landing  gear.  If  the  loads  are  symmetric  with 
respect  to  the  z  =  0  plane  (fig.  5a),  the  loading  can  be  represented  by  a  Fou¬ 
rier  cosine  series.  If  for  some  reason  such  as  the  particular  location  of  a 
discontinuity  in  the  pavement  structure  the  z-axis  must  be  taken  in  a  different 
direction  (fig.  5b) ,  the  same  loading  as  in  figure  5a  will  be  represented  by  a 
complete  Fourier  series  consisting  of  both  cosine  and  sine  terms.  It  is  not 
necessary  to  express  the  arbitrary  gear  loading  as  a  Fourier  series  when  the 
general  three-dimensional  solid  analysis  program  (ref.  9)  is  used;  however,  the 
advantage  in  using  the  prismatic  solid  analysis  technique  in  the  AFPAV  code  is 
that  it  is  simpler  and  less  expensive  than  the  general  three-dimensional  solid 
analysis  for  solving  a  majority  of  practical  airfield  pavement  problems  (refs. 
1,2,4,10). 


The  relative  dimensions  of  the  width,  2b,  and  the  depth,  2d,  with  respect  to 
the  length  of  the  prismatic  solid,  21,  determine  the  appropriate  analytical 
method  (e.g.,  if  b/l  «  1  and  d/Ji  «  1,  conventional  beam  theory  is  adequate). 
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Figure  4.  Pavement  Structure  Idealized 
as  Prismatic  Solid 


Figure  6  shows  an  idealization  of  a  typical  pavement  structure  loaded  by 
the  24  wheels  of  the  landing  gear  and  the  4  wheels  of  the  nose  gear  of  a  C-5A 
aircraft.  Figure  7  shows  the  finite  element  grid  pattern  of  this  idealized 
pavement  structure.  In  this  idealization,  it  is  assumed  that  the  cross  section 
of  the  prismatic  solid  is  as  shown  in  figure  7  with  the  width  equal  to  2B.  The 
period  of  the  applied  loading  is  equal  to  2L  ns  shown  in  figure  6.  (The  numer¬ 
ical  values  of  L  and  B  were  chosen  as  described  in  section  III.)  All  the  ac¬ 
tual  infinite  dimensions  of  the  pavement  system  were  assumed  to  be  finite  as 
shown  by  the  dashed  boundary  lines  as  required  by  this  analysis  technique.  Al¬ 
though  all  24  wheels  of  the  main  landing  gear  are  included  in  figures  6  and  7, 
it  is  not  always  necessary  to  consider  all  the  wheels  in  the  analysis.  In  most 
instances,  only  one  group  of  the  closely  spaced  wheels  need  be  included  (fig. 

8)  since  the  interaction  between  the  different  groups  of  loads  is  usually  insig¬ 
nificant.  Furthermore,  the  direction  of  the  longitudinal  axis  can  be  altered 
as  required  by  the  circumstances  of  the  problem  as  in  the  two  cases  shown  in 
figure  8.  For  instance,  case  2  would  be  the  required  idealization  if  a  cul¬ 
vert  running  under  the  runway  were  to  be  analyzed,  since  the  z-axis  should  coin¬ 
cide  with  the  longitudinal  axis  of  the  culvert. 
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(a)  Loading  Symmetric  About  z  =  0  Plane 


Culvert 


(b)  Loading  Unsymmetric  About  z  =  0  Plane 


Figure  5.  Twelve-Wheel  Assembly  of  C-5A  Landing  Gear 


Main  Landing 
Gear  (24  Wheels), 
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Subbase  Course 


Subgrade 


Idealization  of  Pavement  Structure  Loaded  by  C-5A  Aircraft 


Figure  7.  Finite  Element  Grid  Pattern  of 
Idealized  Pavement  Structure 


11 


2 .  TFEORY  OF  AFPAV  CODE 


Crawford  at  NCEL  developed  the  AFPAV  code  for  AFWL  as  a  model  for  analyz¬ 
ing  airfield  pavements  by  adding  a  preprocessor  (AFPRE)  and  a  postprocessor 
(AFPOST)  to  the  prismatic  solid  analysis  program  (AFPAV)  which  is  the  central 
or  main  part  of  the  code.  The  structural  analysis  of  the  problem  is  performed 
by  AFPAV,  and  the  theory  described  in  the  following  paragraphs  pertains  to  that 
part  only.  AFPRE  simplifies  the  enormous  amount  of  input  data  pertaining  to 
the  multi -wheel  gear  loads  and  the  pavement  structure  properties;  AFPOST  prints 
the  output  data  or  produces  a  graphic  display  of  the  output  either  in  a  movie 
or  in  microfilm  plots. 

As  the  term  finite  element  analysis  indicates,  the  layered  pavement  struc¬ 
ture  with  infinite  dimensions  is  first  reduced  (idealized)  to  a  structure  with 
finite  dimensions  as  shown  in  figure  6  and  then  subdivided  into  a  system  of 

discrete,  but  interconnected  substructures  as  shown  in  figure  7.  Each  rcctan- 

* 

gle/quadrilateral  in  figure  7  represents  a  prismatic  element  (substructure) 
with  a  length  of  2L  (fig.  6).  The  corners  of  these  elements  represent  the  nod- 
aL  lines  and  are  simply  called  nodes. 

Structural  analysis  of  the  complicated  layered  pavement  system  is  performed 
by  assuming  that  the  unknowns  in  the  problem  are  the  displacements  at  the  nodes. 
The  applied  loading  is  broken  up  into  concentrated  forces  acting  at  the  nodes. 

By  applying  the  minimum  potential  energy  theorem  and  using  appropriate  kinemat¬ 
ic  assumptions  for  the  distribution  of  displacements  within  each  element  as  des¬ 
cribed  below,  equilibrium  equations  for  nodal  forces  are  developed  in  terms  of 
the  unknown  nodal  displacements.  Solution  of  this  set  of  simultaneous  equations 
gives  the  nodal  displacements  from  which  the  strains  and  stresses  in  the  differ¬ 
ent  elements  of  the  structure  are  determined. 

In  the  first  version  of  the  code,  each  quadrilateral  element  (fig.  9)  was 
divided  into  four  triangular  elements,  and  a  linear  distribution  of  displacement 
field  was  assumed  within  each  triangular  element.  After  forming  the  stiffness 
matrices  of  the  triangular  elements,  the  quadrilateral  element  stiffness  was 


Henceforth,  in  this  report,  the  use  of  the  adjective  prismatic  to  describe  the 
element  is  dispensed  with  for  convenience;  it  must,  however,  be  remembered  that 
in  the  AFPAV  code  the  author  is  always  referring  to  three-dimensional  prismatic 
elements  of  the  prismatic  solid  and  not  to  two-dimensional  plane  elements. 
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Figure  9.  Prismatic  Quadrilateral  Element 
Composed  of  Triangular  Elements 


formed  by  condensing  out  (i.e. ,  by  eliminating  the  unknown  displacements  of) 
the  interior  node. 

Since  the  nodal  displacements  due  to  arbitrary  gear  wheel  loads  are  func¬ 
tions  of  x,  y,  and  z,  the  displacements  are  expanded  into  Fourier  series  having 
a  period  of  2L,  as  given  in  the  following  equations: 
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where  u,  v,  and  w  represent  displacements  in  the  x,  y,  and  z  directions,  re¬ 
spectively  (ref.  2).  The  Fourier  coefficients  u  ,  u  .  etc.,  arc  functions 
of  x  and  y  only,  and  the  displacement  variation  along  the  z-axis  is  accounted 
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for  by  the  Fourier  series  expansion  in  terms  of  z.  If  the  loading  is  symmet¬ 
rical  about  the  x-axis  (i.e.,  z  =  0  plane),  the  second  terms  on  the  right  side 
of  eqs.  (1)  are  omitted  since  they  are  required  for  unsymmetrical  loading  only. 

The  following  form  of  the  Fourier  coefficients  can  be  assumed  to  give  a 
linear  variation  of  displacements  within  the  triangular  elements  in  the  x-y 
plane : 

u  (x,y)  *  a,  +  a,  x  +  a  y  (2a) 

xn1,  '  in  2n  3n7  y  J 

u  (x,y)  -  a  +  a  x  +  a  y  (2b) 

ynv  i*n  5n  6n7  v  ' 

u  (x,y)  =  a  +  a  x  +  a  y  (2c) 

znv  ?n  sn  an7  v  J 

The  as  arc  called  generalized  coordinates.  The  linear  displacement  field  in 
eqs.  (2)  ensures  a  linear  variation  of  displacements  along  the  edges  of  the  tri¬ 
angles  thus  forcing  continuity  between  elements  since  edges  initially  straight 
will  remain  straight  in  the  deformed  position.  Furthermore,  displacement  con¬ 
tinuity  throughout  the  structure  is  ensured  by  having  common  displacements  at 
the  two  node  points  along  each  element  edge. 

From  the  assumed  displacement  functions,  the  element  strains  are  derived 
by  using  the  appropriate  strain/displacement  relationships  of  elasticity  theory. 
Strains  must  also  be  expanded  in  Fourier  series  since  displacements  are  express¬ 
ed  in  this  manner.  By  using  the  appropriate  constitutive  equations,  the  element 
stresses  can  be  derived  from  the  element  strains.  Stresses,  like  strains,  arc 
expanded  in  Fourier  series. 

The  equilibrium  equations  of  the  structure  are  derived  by  applying  the 
theory  of  minimum  potential  energy.  In  this  process,  when  the  total  strain  en¬ 
ergy  is  computed  by  integrating  over  the  volume  of  the  body,  the  orthogonality 
conditions  of  the  trigonometric  functions  help  to  uncouple  eqs.  (1)  (ref.  2, 
appendix  I).  Thus  a  set  of  governing  equilibrium  equations  in  the  following 
form  is  obtained  for  each  Fourier  term: 

[K  ]{U  }  =  (F  }  (3) 

where  [Kn]  is  the  structure  stiffness  matrix,  {FnJ  is  the  global  nodal  force 
vector,  and  (U  }  is  the  nodal  displacement  vector,  for  the  nth  harmonic.  The 
displacement  vector  in  cq.  (3)  is  a  function  of  x  and  y  only.  The  structure 
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stiffness  matrix  is  obtained  by  summing  the  quadrilateral  element  stiffness 
matrices;  the  quadri lateral  element  stiffness  matrix  itself  is  obtained  from 
the  four  triangle  stiffness  matrices. 

The  discrete  wheel  loads  of  the  landing  gear  of  an  aircraft  are  expressed 
in  a  Fourier  series  as  functions  of  the  z -coordinate  in  the  same  manner  as  the 
displacement  functions.  These  Fourier  loads  are  broken  down  into  an  equiva¬ 
lent  set  of  concentrated  loads  acting  at  the  nodal  points.  In  a  general  case, 
both  cosine  and  sine  Fourier  terms  are  included  in  the  Fourier  series  expansion 


where  f  and  f  are  the  Fourier  coefficients  of  the  applied  vertical  loading 
yn  ym  ' 1  n 

function.  In  eq.  (4),  the  subscript  y  refers  to  the  vertical  component  of  the 
applied  load.  If  loads  in  the  horizontal  direction  arc  to  be  considered,  simi¬ 
lar  expressions  can  be  written.  If  the  wheel  loads  are  symmetric  about  the  x- 
axis  (i.e.,  z  =  0  plane),  eq.  (4)  can  be  reduced  to 
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The  choice  of  the  coordinate  axes  shown  in  figure  6  and  in  case  1  of  figure  8 
requires  the  use  of  eq.  (4);  case  2  in  figure  8  needs  only  eq  (5). 


For  each  harmonic  of  the  loading  function,  there  are  a  corresponding  har¬ 
monic  of  the  displacement  function  and  an  associated  stiffness  matrix.  The 
equilibrium  equations  are  solved  for  each  harmonic  to  obtain  the  Fourier  coef¬ 
ficients  of  the  corres|Mjnding  harmonic  of  the  displacement  function.  The  nodal 
displacements  arc  finally  computed  by  substituting  these  coefficients  in  eqs. 
(1).  The  strains  arc  derived  from  the  displacements.  The  linear  displacement 
functions  assumed  in  eqs.  (2)  give  constant  strains  in  a  given  clement.  1110 
stresses  arc  then  computed  from  the  strains.  Since  the  stresses  arc  constant 
over  a  given  (constant-strain)  clement,  great  care  must  be  exercised  in  ideal¬ 
izing  the  pavement  structure  by  having  smaller  elements  in  areas  of  high-strain 
gradients  and  having  as  many  elements  in  the  whole  structure  as  computer  capac¬ 
ity  and  computer  cost  will  permit.  However,  because  of  the  limitations  in  com¬ 
puter  capacity,  it  is  not  possible  to  increase  the  number  of  elements  too  much; 
therefore,  it  is  necessary  to  replace  the  constant -strain  elements  by  higher- 
order  elements  where  greater  accuracy  is  desired  (ref.  11).  Details  of  the 
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highei  rdcr  elements  and  a  comparison  of  the  pavement  responses  predicted  by 
constant-strain  clement  formulation  and  higher-order  element  formulation  will 
be  given  in  a  subsequent  report. 

3.  STRUCTURE  OF  AFPAV  CODE 

The  AFPAV  code  contains  three  parts  each  of  which  performs  a  significant 
function:  AFPRE  (preprocessor)  generates  the  finite  element  mesh  and  deter¬ 
mines  the  Fourier  coefficients  of  the  applied  wheel  loading  for  input  into  the 
main  program;  AFPAV  (main  program)  uses  the  linear  elasticity  theory  to  predict 
the  response  of  layered  pavement  systems  subjected  to  the  static  loadings  of 
MW1IUL  of  aircraft  such  as  the  C-5A  and  the  Boeing  747;  AFPOST  (postprocessor) 
is  capable  of  either  producing  a  graphic  display  of  the  stresses  and  displace¬ 
ments  or  of  printing  out  the  stresses,  strains,  and  displacements  in  the  entire 
cross  section  of  the  structure  at  any  desired  location  in  the  pavement. 

a.  AFPRE  Program 

The  main  program  of  the  AFPAV  code  considers  the  effect  of  multiple 
wheel  loads  without  converting  them  into  an  equivalent  single-wheel  load  or 
without  using  the  principle  of  superposition  as  in  the  BISTRO  code  (ref.  5). 
Furthermore,  the  Air  Force  has  developed  this  code  as  an  analytical  model  to 
represent  the  special  three-dimensional  nature  of  the  pavement  structure  (i.e., 
prismatic  solid)  and  to  consider  the  actual  gear  configuration  of  any  type  of 
aircraft.  Because  of  this  versatility,  a  large  amount  of  manual  input  will  be 
required  to  define  the  pavement  structure  and  the  Fourier  coefficients  of  ap¬ 
plied  loading  if  only  the  original  prismatic  solid  finite  element  analysis  pro¬ 
gram  developed  by  Herrmann  is  used.  Therefore,  AFPRE  was  developed  to  automate 
the  input  of  structure  dimensions  and  the  calculation  of  the  Fourier  coeffi¬ 
cients  of  the  applied  loading. 

AFPRE  contains  several  options  and  special  features.  There  arc  two 
distinct  segments  in  the  preprocessor  (viz.,  single-wheel  and  multi-wheel  load 
analysis).  In  the  single-wheel  analysis,  only  one-half  of  the  layered  pavement 
system  (fig.  10)  need  be  considered.  The  main  restriction  in  this  analysis  is 
that  the  thickness  of  each  layer  of  the  pavement  system  be  uniform.  For  case 
in  expressing  the  wheel  loads  in  Fourier  series,  the  tire  contact  area  is  as¬ 
sumed  to  be  rectangular.  The  multi-wheel  load  option  in  AFPRE  is  more  versa¬ 
tile  than  the  single-wheel  option.  There  arc  two  cases  in  the  multi-wheel  load 


option,  depending  on  the  choice  of  the  z-axis  with  respect  to  the  gear  wheel 

configuration  in  a  given  problem.  They  are  symmetric  and  unsymmetric  wheel 

configuration.  Figure  11  shows  a  6-wheel  assembly  symmetric  with  respect  to 
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the  longitudinal  z-axis.  Because  of  symmetry,  it  is  sufficient  to  consider 
only  one-half  of  the  pavement  system  loaded  by  one-half  of  the  wheel  assembly. 
This  feature  is  similar  to  the  single -wheel  load  option  in  AFPRE  and  reduces 
the  number  of  node  points  and  elements.  Of  course,  this  assumes  symmetry  of 
not  only  the  loads  but  also  of  the  pavement  layers.  The  same  6-wheel  assembly 
is  shown  again  in  figure  12,  but  the  orientation  of  the  z-axis  is  different 
from  that  in  figure  11  because  of  the  inclusion  of  the  culvert  in  the  pavement 
structure.  Therefore,  the  entire  wheel  assembly  must  be  included  in  the  anal¬ 
ysis. 

There  is  a  significant  difference  between  the  single-wheel  option  and 
the  multi -wheel  option  in  respect  to  the  capability  for  automatic  mesh  genera¬ 
tion.  In  the  single-wheel  option,  the  finite  element  mesh  data  for  the  entire 
structure  is  automatically  generated;  however,  to  achieve  this  capability,  a 
sacrifice  must  be  made  in  that  the  pavement  layers  must  be  horizontal  and  of 
uniform  thickness.  If  nonuniform  thicknesses  are  encountered  in  a  problem,  the 
multi -wheel  load  option  can  be  used  to  input  the  structure  data  even  for  a  sin¬ 
gle  wheel  load.  In  the  multi -wheel  option,  the  following  three  distinct  re¬ 
gions  are  assigned  to  cover  the  different  layers  of  the  pavement  systems: 

(1 j  Surface  region 

(2)  Base  region 

(T)  Subgrade  region 

liven  though  the  names  of  the  different  regions  appear  to  be  the  same 
as  those  commonly  used  by  pavement  engineers,  they  do  not  necessarily  refer  to 
the  same  thing  in  the  AFPAV  code.  For  example,  in  figure  11  there  are  four  dif¬ 
ferent  pavement  material  layers,  viz.,  surface  course,  base  course,  subbasc, 
and  subgrade.  In  each  of  these  material  layers,  there  arc  more  than  one  clement 
layer.  In  the  AFPAV  code,  however,  the  first  of  the  clement  layers  of  the  sur¬ 
face  course  is  called  the  surface  region;  the  remaining  clement  layers  of  the 


However,  the  loading  is  not  symmetric  about  the  z  =  0  plane  in  this  case  and, 
therefore,  eq.  (4)  containing  both  cosine  and  sine  Fourier  terms  would  be  used 
to  determine  the  Fourier  coefficients  of  the  loading. 
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surface  course  and  all  the  clement  layers  of  the  base  course  and  subbase  are 
designated  as  the  hue  region ;  and  the  subgrade,  which  is  always  assumed  to  be 
of  one  homogeneous  material  in  this  code  for  the  purpose  of  automation  of  ele¬ 
ment  layers,  is  called  the  uubgrade  region.  This  division  of  the  different 
pavement  layers  into  three  regions  has  been  arbitrarily  made  to  achieve  very 
desirable  objectives.  While  inputting  the  information  in  AFPRE  concerning  dif¬ 
ferent  layers  for  generating  the  finite  element  mesh,  it  is  necessary  to  give 
the  coordinates  of  the  nodes  and  the  node  numbers  of  the  elements  in  the  sur¬ 
face  region  only.  Then,  if  the  number  of  element  layers  required  in  each  oT  the 
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component  material  layers  of  the  base  region  and  the  thicknesses  of  the  materi¬ 
al  layers  arc  specified,  AFPRE  automatically  calculates  all  the  information  con¬ 
cerning  the  elements  and  the  nodes  in  the  entire  structure.  The  sides  of  the 
surface  and  the  base  regions  arc  assumed  to  be  vertical ;  the  subgradc  region 
has  sloping  sides. 

In  some  special  cases  the  surface  region,  like  the  base  region,  can 
have  more  than  one  material.  (This  can  be  seen  in  figure  12  where  the  surface 
course,  the  base  course,  the  culvert,  and  even  a  part  of  the  subgrade  arc  in¬ 
cluded  in  the  surface  region.)  This,  however,  increases  the  amount  of  input 
data  to  describe  the  complicated  structure  to  AFPRE. 

In  the  present  version  of  AFPRE,  the  nodes  arc  first  numbered  auto¬ 
matically  in  a  horizontal  direction.  Then,  if  desired,  the  nodes  can  be  renum¬ 
bered  in  the  vertical  direction  as  shown  in  figure  13.  This  results  in  a  re¬ 
duction  of  the  computer  time  required  for  the  solution  of  the  set  of  simultan¬ 
eous  equations  represented  by  cq.  f3)  because  of  the  reduction  in  the  bandwidth 
of  the  matrix  equation.  I igure  13  shows  the  finite  element  idealization  of  a 
typical  [)avement  system  with  14  element  layers  (and  3  or  4  material  layers  as 
the  case  may  be).  Mien  the  numbering  of  the  40  columns  of  nodes  is  done  hori¬ 
zontally,  the  bandwidth  in  this  three-dimensional  problem  is  126.  However,  by 
renumbering  the  nodes  in  a  vertical  direction  as  shown  in  figure  13,  the  band¬ 
width  is  reduced  to  51.  Since  the  number  of  arithmetic  operations  involved  in 


The  number  of  element  layers  in  the  subgradc  region  is  automatically  calcu¬ 
lated;  however,  the  user  must  specify  certain  information  to  guide  the  program 
in  determining  the  depths  of  the  clement  layers. 
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Figure  13.  Renumbering  of  Nodes  of  Idealized  Pavement  Structure 


solving  a  set  of  N  simultaneous  equations  with  a  bandwidth,  M,  is  proportional 
to  ,NM2,  the  renumbering  reduces  the  number  of  operations  to  about  16  percent  of 
that  required  without  renumbering.  (If  the  number  of  element  layers  were  20 
instead  of  14,  the  reduction  in  the  number  of  operations  would  be  about  30  per¬ 
cent.  The  saving  due  to  renumbering  is  obviously  a  function  of  the  number  of 
horizontal  layers  for  a  given  number  of  vertical  columns.)  This  results  in 
considerable  savings  in  computer  time,  particularly  when  solving  a  multi-layered 
pavement  system  loaded  by  multi -wheel -geared  aircraft  such  as  the  C-5A.  At 
present,  this  renumbering  scheme  requires  the  running  of  AFPRE  twice;  however, 
AFPRE  could  be  modified  to  eliminate  this  double  run.  (AFPRE  takes  a  very  in¬ 
significant  part  of  the  total  computer  time  compared  with  that  taken  by  the  com¬ 
plete  AFPAV  code.) 

b.  AFPAV  Main  Program 

The  main  part  of  the  AFPAV  code  takes  the  input  data  concerning  the 
aircraft  loading  and  the  pavement  structure  from  APPRE  and  performs  the  finite 
element  structural  analysis  computations.  In  the  present  version  of  this  code, 
this  program  performs  only  static  analysis  of  linear  elastic  problems;  however, 
this  part  of  the  code  is  being  modified  to  consider  nonlinear  stress/strain  re¬ 
lationships  of  certain  pavement  layer  materials. 

There  are  two  options  for  feeding  the  input  data  to  the  main  program: 
(1)  using  AFPRE  for  generating  the  finite  clement  mesh  data  and  determining  the 
Fourier  coefficients  of  the  applied  loading;  and  (2)  determining  the  input  data 
by  separate  calculations  and  feeding  them  to  the  main  program  manually  by  cards. 
Since  the  purpose  of  adding  AFPRE  to  the  original  prismatic  solid  analysis  pro¬ 
gram  was  to  eliminate  the  need  for  manual  input  of  an  enormous  amount  of  data 
relating  to  the  finite  element  mesh  and  multiple  loads,  it  is  not  now  necessary 
to  use  the  second  input  option.  Furthermore,  AFPRE  does  not  take  much  computer 
tunc. 

There  are  two  options  for  obtaining  the  output  data:  (1)  using  AFKIST; 
and  (2)  using  the  main  program  only.  Using  the  main  program  to  print  the  out¬ 
put,  one  can  specify  one  or  more  nodes  or  elements  at  one  or  more  cross  sections 
of  the  pavement  system  and  the  following  information  pertaining  to  any  specified 
cross  section  of  the  pavement  system  can  be  obtained: 

(1)  The  displacements  in  three  directions  at  any  speci¬ 
fied  node  of  the  finite  clement  mesh 
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(2)  The  three  normal  strains,  the  three  shear  strains, 
and  the  three  normal  stresses  in  any  specified 
element 

Using  AFPOST,  however,  one  can  limit  the  number  of  cross  sections  but  not  the 
number  of  elements  and  nodes  in  a  given  cross  section.  Once  the  cross  sections 
are  specified,  the  results  relating  to  all  the  nodes  and  the  elements  in  those 
cross  sections  are  printed.  When  using  AFPOST,  the  core  storage  required  in¬ 
creases  depending  on  the  number  of  cross  sections  at  which  results  are  desired. 

c .  AFPOST  Program 

To  eliminate  the  need  to  examine  the  enormous  amount  of  printed  out¬ 
put  generated  either  by  AFIXXST  or  by  the  main  program  itself,  certain  plot  cap¬ 
abilities  were  developed  to  do  the  following  for  any  specified  cross  section: 

(1)  Plot  contours  of  stresses  and  displacements  in  the 
x-y  plane 

(2)  Plot  stresses  and  displacements  as  surfaces  over 
the  x-y  plane  providing  a  three-dimensional  view 
of  the  distribution  of  stresses,  etc.,  in  a  cross 
section 

These  plots  can  be  obtained  for  the  entire  cross  section  of  the  pavement  system 
or  for  only  a  part  of  the  cross  section  as  desired  by  the  user. 

In  addition  to  the  above  capabilities,  AFPOST  can  plot  the  finite  elc 
meat  mesh  input  data  (such  as  the  material  description,  mesh  geometry,  node  num 
hers,  etc.j  so  that  the  user  can  check  the  accuracy  of  the  input,  hxamplcs  of 
most  of  the  plots  that  can  be  obtained  using  AFPOST  arc  given  in  reference  2. 

AFPOST  can  also  help  the  user  make  a  movie  showing  all  the  plots  men¬ 
tioned  above.  This  movie  does  not  show  the  variation  of  stresses,  etc.,  as  a 
function  of  time,  since  only  static  analysis  can  presently  be  performed.  How¬ 
ever,  it  does  afford  an  excellent  graphic  display  of  the  changes  in  stresses 
and  displacements  as  a  function  of  space.  In  the  case  of  a  C-5A  aircraft,  this 
graphic  display  could  be  assumed  to  represent,  without  any  significant  error, 
the  changes  in  stresses,  etc.,  at  a  given  cross  section  of  the  pavement  system 
as  the  different  rows  of  the  24-whcel  assembly  of  the  landing  gear  moved  very 
slowly  over  that  cross  section.  Such  a  movie  has  been  made,  and  prints  of  the 
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frames  will  be  included  in  a  subsequent  report  which  will  contain  documenta¬ 
tion  of  the  program  and  sample  problems. 

In  addition  to  the  graphic  displays  mentioned,  AFI’OST  can  print  out 
in  tabular  form  all  the  stresses  and  displacements  in  all  of  the  elements  at 
desired  cross  sections  of  the  pavement  structure  as  stated  earlier.  However, 
while  it  is  possible  to  print  out  displacements,  strains,  and  stresses  pertain¬ 
ing  to  only  a  select  number  of  nodes  or  elements  using  the  main  program,  it  is 
not  possible  to  achieve  this  using  the  present  version  of  AFPOST.  Experience 
has  shown  that  it  is  not  advisable  to  use  AFPOST  if  the  pavement  response  at 
only  a  few  points  is  desired.  For  graphic  displays  of  contours  of  stresses, 
etc . ,  however ,  AFPOST  is  very  useful . 

The  preceding  description  of  the  capabilities  of  AFPOST  covers  only 
those  features  originally  included  in  the  code  developed  by  NCEL.  While  analyz 
ing  practical  pavement  problems,  it  was  realized  that  it  would  be  useful  if  lin 
car  plots  of  stresses  and  displacements  as  functions  of  depth  could  also  be 
made.  Therefore,  a  subroutine  for  obtaining  such  plots  was  written  at  CERF  by 
Ann  Stough.  This  subroutine  can  be  used  in  conjunction  with  the  main  program 
of  the  AFPAV  code.  Similarly  it  would  be  useful  to  have  the  capability  to  plot 
the  surface  deflection  basins  of  pavement  systems  loaded  by  multi -wheeled  air¬ 
craft  gears  using  the  output  data  given  by  the  main  program.  A  subroutine  to 
obtain  such  a  plot  was  also  prepared  at  CERF.  In  order  not  to  make  AFPOST  a 
very  unwieldy  part  of  the  AFPAV  code,  these  two  subroutines  arc  kept  separate 
and  can  be  used  in  conjunction  with  the  main  program  whenever  required. 

4.  SPECIAL  FEATURES  AND  LIMITATIONS 

A  special  feature  of  the  AFPAV  code  is  that  the  culvert  shown  in  figure  12 
is  considered  a  part  of  the  pavement  structure  in  the  prismatic  solid  analysis; 
this  cannot  be  done  with  other  codes  used  for  pavement  analysis.  Similarly, 
edge  loading,  discontinuities  in  pavement  systems  such  as  longitudinal  joints, 
and  nonuniform  thickness  of  certain  layers  (e.g.,  thickened  edges  and  keel  sec¬ 
tions  in  Portland  cement  concrete  slabs)  can  be  considered  by  the  AFPAV  code 
and  not  by  other  codes  such  as  BISTRO.  For  a  detailed  comparison  of  the  AFPAV 
code  with  other  codes,  see  reference  4. 

Although  the  AFPAV  code  is  a  good  tool  for  three-dimensional  analysis  of 
pavement  systems  idealized  as  periodically  loaded  prismatic  structures,  a 
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general  three-dimensional  analysis  of  pavement  structures  cannot  be  performed 
with  this  code.  Specifically,  the  presence  of  intersecting  joints  in  Portland 
cement  concrete  slabs  cannot  be  accounted  for,  since  the  main  assumption  of  a 
prismatic  solid  is  that  the  cross  section  of  the  solid  does  not  vary  in  the  lon¬ 
gitudinal  direction.  In  the  present  version  of  the  AFPAV  code,  only  linear 
clastic,  static  analysis  can  be  performed.  However,  this  code  is  currently  be¬ 
ing  modified  to  consider  nonlinear  stress/strain  relationships  of  certain  pave¬ 
ment  layer  materials. 
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SKCTION  III 


SIGNIFICANCE  OF  SYSTEM  IDEALIZATION  PARAMETERS 

Both  the  pavement  structure  and  the  discrete  loads  of  the  multi -wheeled 
aircraft  gear  must  be  properly  idealized  in  order  to  analyze  a  layered  air¬ 
field  pavement  problem.  The  structure  idealization  of  pavement  systems  for 
finite  element  analysis  has  been  discussed  in  detail  elsewhere  (refs.  1,2,4, 
10,12,13).  Briefly,  the  layered  half-space  is  idealized  as  a  finite  structure 
with  boundaries  consisting  of  rollers  and  pin  supports  as  shown  in  figures  10 
through  13,  and  the  discrete  loads  are  expressed  in  a  Fourier  series.  The  in¬ 
fluence  of  these  boundary  conditions  and  the  significance  of  the  loading  func¬ 
tion  parameters  are  discussed  in  this  section.  The  two  parameters  pertaining 
to  the  Fourier  series  expansion  of  the  loading  function  are  (1)  the  number  of 
Fourier  terms,  N,  and  (2)  the  characteristic  length,  or  the  period,  2L,  as  in 
eq.  (1).  Linked  with  these  two  parameters  is  the  fineness  of  the  finite  ele¬ 
ment  mesh  used  to  represent  the  given  pavement  cross  section.  The  interaction 
of  these  three  parameters  has  been  investigated  by  Herrmann  whose  findings  on 
this  subject  are  discussed  by  Crawford  in  some  detail  in  reference  2.  In  this 
connection,  it  is  worthwhile  to  make  the  following  general  observations  based 
on  experience  gained  in  using  the  code. 

1.  ELEMENT  SIZE 

Crawford  (ref.  2)  has  stated  that,  when  the  decreasing  importance  of  the 
higher-order  terms  to  the  total  solution  is  accounted  for,  about  the  same  fine¬ 
ness  of  grid  is  needed  for  all  values  of  N.  Therefore,  it  is  sufficient  to 
draw  up  a  finite  element  grid  with  element  sizes  under,  and  adjacent  to,  the 
loads  as  small  as  the  computer  capacity  will  permit  and  to  use  the  same  grid 
for  all  values  of  N. 

Figure  13  shows  a  finite  element  mesh  used  to  determine  the  response  of  a 
pavement  system  loaded  by  a  C-5A  aircraft.  The  distance  of  the  vertical  roll¬ 
ers  from  the  edge  of  the  outer  wheels,  x,  should  be  as  great  as  possible  so  as 
to  include  as  much  of  the  pavement  cross  section  in  the  idealized  structure  as 
is  convenient.  This  distance  should  be  at  least  7  to  10  times  the  width  of  the 
wheel,  and  preferably  greater  in  the  case  of  rigid  pavement  systems.  It  may  be 
difficult  to  achieve  this  for  multiple  wheels  because  of  the  increase  in  the 
number  of  nodes  that  will  result  from  an  excessively  large  value  of  x.  If  in 


28 


figure  13  the  top  element  layers  are  taken  as  1  in.  thick  in  a  3-in. -thick 
material  such  as  an  asphalt  cement  concrete  layer,  it  would  he  desirable  to 
restrict  the  maximum  width  of  an  element  in  this  clement  layer  at  either  edge 
of  the  structure,  to  about  twenty  times  the  thickness  (i.c.,  20  in.)  to  obtain 
a  proper  definition  of  stresses  in  that  element.  If  this  restriction  is  ad 
hcred  to  and  if  the  distance,  x,  chosen  is  200  in.  or  more,  the  number  of  nodes 
in  a  horizontal  layer  will  increase  to  considerably  more  than  that  shown  in 
figure  13.  The  increase  in  the  number  of  nodes  in  a  horizontal  layer  will  re¬ 
strict  the  number  of  element  layers  since  there  is  a  limit  on  the  total  number 
of  nodes  that  the  program  can  handle.  Since  the  computer  time  depends  on  the 
number  of  nodes  in  the  grid,  it  is  desirable  to  minimize  the  number  of  nodes  in 
a  given  problem.  Therefore,  care  should  be  exercised  in  choosing  the  finite 
element  grid  to  idealize  a  given  pavement  system  loaded  by  a  particular  air¬ 
craft  . 

2.  MJMRIiR  Of  FOUR  I  HR  TF.RMS 

The  number  of  Fourier  terms,  N,  required  to  express  a  loading  function 
with  reasonable  accuracy  must  be  determined  in  relation  to  the  problem  to  bo 
solved.  In  reference  4  the  author  shows  that  the  clastic  (recoverable!  sur¬ 
face  deflections  are  predicted  fairly  accurately  in  a  particular  problem  with 
as  little  as  six  Fourier  terms.  Although  the  subgrade  stresses  were  predicted 
as  accurately  with  6  Fourier  terms  as  with  31,  the  stresses  in  the  surface  lay¬ 
er  were  not  predicted  properly  with  a  small  number  of  Fourier  terms.  The  rea¬ 
son  for  this  is  clear  from  figure  14  which  shows  the  Fourier  series  representa¬ 
tion  of  a  line  of  four  wheel  loads  of  a  C-5A  gear  by  9,  17,  and  31  Fourier 
terms.  Figure  15,  which  compares  the  loading  function  represented  by  9,  17, 
and  31  terms,  also  shows  that  a  smaller  number  of  Fourier  terms  docs  not  pre¬ 
dict  the  surface  stresses  accurately  since  large  tensile  (upward)  loading  is 
given  by  Fourier  expansion  with  a  smaller  N.  The  stresses  in  deep  layers,  how¬ 
ever,  will  be  accurately  predicted  as  can  hr*  explained  by  the  Saint-Vcnant 's 
principle.  (A  unit  value  is  assumed  in  figures  14  and  15  for  contact  pres¬ 
sures,  and  a  constant  period  is  used  for  the  three  different  Fourier  terms.) 

3.  PF.RIOD  OF  LOADING  FUNCTION 

The  period  of  the  Fourier  series  representing  the  discrete  wheel  loads  af¬ 
fects  the  pavement  response  just  as  the  number  of  Fourier  terms  does.  Figure 
16  shows  the  effect  of  the  period  on  the  loading  function  of  a  line  of  four 
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Figure  15.  Effect  of  Number  of  Fourier  Terms  on  Loading  Function 
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Figure  lb.  Effect  of  Period  on  Fourier  Series  Representation 
of  Four  Wheel  Loads  of  C-141A  Landing  Gear 


*  ijMjiivujiiiuumnuu 


wheels  of  a  0-141A  aircraft.  When  the  half-period  is  210  in.,  which  is  equal 
to  the  tread,  i.c.,  the  distance  between  the  centroids  of  the  two  sets  of  four 
wheels  ol  the  main  landing  gear  (ref.  14),  the  maximum  value  of  the  loading 
function  (with  N  *  20)  is  0.986  which  is  almost  the  assumed  contact  pressure 
of  1.0.  However,  when  the  half-period  is  increased  to  300  in.  and  400  in., 
the  corresponding  maximum  values  of  the  loading  function  (with  N  =  20)  become 
0.766  and  0.502.  This  adversely  affects  the  prediction  of  surface  layer  stress¬ 
es.  Figure  17  shows  a  comparison  of  the  loading  functions  for  the  wheel  loads 
of  the  C-141A  aircraft  for  the  three  values  of  2L.  As  L  increases,  an  increas¬ 
ing  number  of  Fourier  terms  are  required  to  obtain  a  better  approximation  of 
the  loading  function;  but,  since  this  increases  computer  time,  judgment  based 
on  experience  in  using  AFPRF.  should  be  exercised  in  choosing  an  appropriate 
characteristic  length  for  a  given  problem. 

Figures  18  and  19  show,  respectively,  the  effect  of  the  half-period  on  the 
elastic  surface  deflection  and  on  the  horizontal  stresses  at  the  bottom  (i.e., 
at  a  depth  of  7.2  in.)  of  an  8-in.  Portland  cement  concrete  slab  loaded  by  the 
eight  wheels  of  the  C-141A  aircraft  landing  gear.  There  is  almost  no  differ¬ 
ence  in  the  maximum  surface  deflection  for  the  various  half -periods,  although 
there  is  an  increase  of  about  12  percent  in  the  maximum  horizontal  stress  (fig. 
19) .  Even  though  this  percentage  is  small ,  it  may  be  important  to  know  the  cor¬ 
rect  maximum  tensile  stress  if  it  approaches  the  modulus  of  rupture. 

The  effects  of  the  period  of  the  loading  function  on  flexible  pavement  re¬ 
sponse  may  differ  quantitatively  from  those  on  the  rigid  pavement  response. 

But  even  in  flexible  pavement  systems ,  only  the  response  of  the  upper  layers  is 
affected  by  the  number  of  Fourier  terms  and  the  period  of  the  loading  function, 
which  is  also  called  the  oharaateviatia  length.  It  is,  therefore,  necessary  to 
exercise  reasonable  judgment  in  selecting  proper  values  of  N  and  L  when  analyz¬ 
ing  a  given  problem,  keeping  in  mind  the  nature  of  the  information  required. 

As  stated  earlier,  if  the  total  surface  deflection  is  all  that  is  required  from 
an  analysis,  any  set  of  reasonable  N  and  L  values  will  serve  the  purpose,  and 
no  time  need  be  wasted  in  representing  the  loading  function  very  accurately. 

4.  BOUNDARY  CONDITIONS 

Another  important  consideration  when  idealizing  a  layered  pavement  struc¬ 
ture  is  the  selection  of  an  appropriate  width  for  the  idealized  structure. 
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Figure  IS.  Effect  of  Period  on  Elastic  Surface  Deflection  of  Rigid  Pavement 


Figure  19.  Effect  of  Period  on  Maximun  Horizontal 
Stress  of  Rigid  Pavement 


1‘his  is  just  as  significant  as  the  selection  of  the  characteristic  length.  It 
would  be  ideal  if  the  entire  width  of  the  airfield  pavement  could  be  used  when 
drawing  up  the  finite  element  mesh;  but,  this  is  impossible  because  of  the  lim¬ 
itation  on  the  number  of  nodes  and  elements  imposed  by  computer  capacity  and 
cost.  The  problem  of  choosing  a  proper  width  for  the  idealized  structure  is, 
in  fact,  a  part  of  the  overall  problem  of  the  boundary  conditions  assumed  at 
the  two  sides  and  at  the  bottom  of  the  idealized  pavement  structure.  As  shown 
in  figures  10  through  13,  the  bottom  supports  are  assumed  to  allow  no  transla¬ 
tional  movement  in  either  the  x-  or  y-dircction,  and  these  supports  could  be  at 
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a  depth  of  fifty  to  sixty  times  the  half-width  of  a  wheel  contact  area.  The 
sides  are  assumed  to  consist  of  roller  supports  which  restrain  movement  perpen¬ 
dicular  to  the  sides.  The  surface  and  base  regions  have  vertical  sides;  the 
subgradc  region  has  sloping  sides  (angle  of  slope  *  30°  to  the  vertical  axis). 
The  width  of  a  flexible  pavement  structure  could  be  obtained  by  assuming  that 
the  vertical  rollers  are  at  a  distance  of  about  20  times  the  half-width  of  the 
outer  wheel  contact  area;  but,  the  width  of  a  rigid  pavement  structure  should 
be  obtained  by  placing  the  vertical  rollers  as  far  away  as  possible  from  the 
outer  wheel,  say  about  30  times  the  half -width  of  the  wheel  contact  area.  In 
this  connection,  it  is  noteworthy  that  the  AFPAV  code  has  an  important  advan¬ 
tage  over  the  WIL67  and  BISTRO  codes  which  perform  axisynmetric  analyses  of 
pavement  systems.  Whereas  BISTRO  and  W1L67  cannot  consider  edge  loading  in  the 
case  of  a  rigid  pavement  system,  the  AFPAV  code  can  analyze  edge  loading  be¬ 
cause  it  treats  the  structure  as  a  prismatic  solid. 

An  eight-wheel  assembly  of  the  main  landing  gear  of  a  0 14 1 A  aircraft  is 
assumed  to  load  a  pavement  with  the  two  sets  of  four  wheels  parallel  to  the 
edge  of  the  pavement  as  shown  in  figure  20.  Figure  21  shows  the  effect  of  the 
location  of  the  load  on  the  normal  stresses  at  some  selected  points  in  an  8- in. 
Portland  cement  concrete  slab  resting  on  a  10-CBR  subgradc.  The  contact  pres¬ 
sure,  p,  is  assumed  to  be  180  psi.  The  boundary  conditions  with  concrete  ex¬ 
tending  up  to  the  vertical  rollers  on  both  sides  (fig.  21a)  represent  an  in¬ 
terior  loading  condition  (i.c.,  one  with  the  aircraft  wheels  located  in  the  in¬ 
terior  part  of  a  wide  airfield  pavement).  Figure  21b  shows  an  edge  loading 

* - 

If  a  rigid  layer  is  known  to  exist  at  a  lesser  depth,  the  actual  thickness  of 
the  subgradc  should  be  used  in  the  finite  clement  idealization  of  the  structure. 


37 


Figure  20.  Location  of  Main  Landing  Gear  of  C-141A 
Aircraft  for  Boundary  Condition  Problem 

condition  (i.c.,  one  with  the  set  of  rear  wheels  coming  to  the  very  edge  of  the 
slab).  The  top  of  the  slab  is  assumed  to  be  at  the  same  level  as  the  adjacent 
ground  which  has  the  same  li-modulus  as  the  subgrade  beneath  the  slab. 

There  is  a  very  significant  difference  in  stress  distribution  between  the 
interior  loading  and  the  edge  loading  conditions.  In  the  former  case,  the  max¬ 
imum  tensile  stress  in  the  longitudinal  direction,  oz,  is  about  500  psi  (point 
B  in  fig.  21a).  In  the  edge  loading  case,  the  maximum  tensile  stress  in  the 
longitudinal  direction,  o  ,  is  about  630  psi  (point  A  in  fig.  21b).  Consider¬ 
ing  the  stresses  at  all  the  points  in  the  slab,  there  is  an  increase  of  about 
26  percent  in  the  maximum  horizontal  tensile  stress.  Considering  only  the 
stresses  at  an  interior  point  such  as  B  under  the  front  set  of  four  wheels, 
there  is  no  significant  difference  in  the  magnitude  of  normal  stresses  (fig. 

21  j;  but,  the  stresses  at  point  A  near  the  edge  of  the  slab  change  dramatical  ly, 
c.g.,  the  value  of  n changes  from  a  compressive  stress  of  6.2  psi  to  a  tensile 
stress  of  630.1  psi.  \gain,  only  the  Al'IV.V  code  can  compute  the  stresses  due 

to  such  edge  loading.  The  BISTRO  code  or  any  other  axisymmctric  code  can  only 
determine  the  stresses  due  to  interior  loading  and,  therefore,  this  code  would 
have  undcrprcdictcd  the  maximum  tensile  stress  caused  by  the  C-141A  aircraft  in 
this  problem. 
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174  in. 


3)  Finite  element  grid  is  not  drawn  to  scale; 
elements  in  the  subgrade  are  not  shown  in 
this  figure. 

(a)  Interior  Loading 

Effect  of  Load  Location  on  Stresses  in  Rigid  Pavement 


Figure  21 ---Concluded 


SECTION  IV 


ANALYSIS  OF  MWHGL  PAVEMENT  TEST  SECTIONS 

A  comprehensive  experimental  research  project  was  jointly  sponsored  by 
the  IJ.S.  Air  Force,  the  II. S.  Army,  and  the  Federal  Aviation  Administration  to 
obtain  data  on  the  behavior  of  different  pavement  systems  subjected  to  multiple 
wheel  heavy  gear  loads.  Flexible  and  rigid  pavement  sections  (fig.  22)  were 
constructed  and  tested  at  IVES  (ref.  15)  to  determine,  among  other  things,  the 
pavement  response  to  (1)  static  loading  and  (2)  trafficking  of  a  prototype  of 
a  12-wheel  assembly  (i.e.,  one  main  landing  gear  of  a  C-5A  aircraft).  The  ex¬ 
perimental  data  obtained  from  this  full-scale  pavement  testing  were  compared 
with  the  theoretical  response  predicted  by  computer  codes. 

Before  the  AFPAV  code  was  fully  developed,  the  BISTRO,  the  WIL67,  and  the 
VISAB3  codes,  along  with  the  main  program  of  the  AFPAV  code  (without  the  prepro¬ 
cessor  and  the  postprocessor)  were  used  in  the  theoretical  analysis  of  these 
test  sections.  The  results  of  this  analysis  were  compared  with  the  experiment¬ 
ally  determined  values  for  the  test  sections  (ref.  4).  After  the  preprocessor 
and  the  postprocessor  were  added  to  the  AFPAV  main  program,  these  test  sections 
were  again  analyzed  using  all  three  parts  of  the  AFPAV  code  as  a  check  on  the 
performance  of  the  completed  AFPAV  code.  While  looking  at  the  comparison  of 
field  observations  and  predicted  data,  the  following  two  important  facts  must 
be  kept  in  mir.vl:  (1)  The  original  version  of  the  AFPAV  code  used  in  this  analy¬ 
sis  was  based  on  a  linear  elastic  model  for  all  of  the  components  of  the  pave¬ 
ment  systems;  and  (2)  the  clastic  properties  of  the  different  layer  materials 
were  estimated  in  the  absence  of  actual  data  obtained  by  procedures  which  would 
give  properties  corresponding  to  the  strain  levels  pertaining  to  the  applied 
loading  (ref.  16). 

1.  FLEXIBLE  PAVEMENT  TEST  SECTION 

As  seen  in  figure  22a,  the  flexible  pavement  test  section  consisted  of  live 
items.  In  all  of  these  items  the  surface  course  was  a  3-in. -thick  asphalt  ce¬ 
ment  concrete  layer  and  the  base  course  was  a  6-in. -thick  crushed  stone  layer; 
the  thicknesses  of  the  sunbasc  and  the  subgrade  of  these  items  varied  (fig. 

22a). 
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Item  4  of  the  flexible  pavement  test  section  was  analyzed  using  the  main 
program  of  the  AFPAV  code  and  the  theoretical  response  was  compared  with  the 
preliminary  experimental  test  data  obtained  at  WES  (ref.  4).  Table  I  presents 
the  layer  properties  that  were  used  for  the  various  layers  of  flexible  item  4 
in  the  analysis  reported  in  reference  4.  It  was  shown  therein  that  the  pave¬ 
ment  response  predicted  by  the  linear  elastic  analysis  did  not  vary  signifi¬ 
cantly  whether  item  4  was  considered  to  consist  of  seven  layers  as  shown  in  fig 
ure  22a  or  to  consist  of  only  four  layers  as  idealized  in  table  I.  Therefore, 
item  3  of  the  flexible  test  section  was  idealized  as  a  four- layered  pavement 
system  and  analyzed.  The  layer  properties  of  item  3  are  shown  in  table  II. 
Since  the  value  of  108,000  psi  for  the  E-modulus  of  the  crushed-stone  base 


Table  I 

LAYER  PROPERTIES  OF  WES  FLEXIBLE  ITB4  4 


Layer 

Thickness , 
in. 

Modulus  of  Elasticity, 
psi 

Poisson's  Ratio 

Surface  Course 

3 

150,000 

0.25 

Base  Course 

6 

108,000 

0.25 

Subbase 

24 

22,500 

0.35 

Subgrade 

00 

4,400 

0.45 

Note:  Data  are  from  reference  4. 


Table  II 

ASSUMED  LAYER  PROPERTIES  OF  WES  FLEXIBLE  ITEM  3 


Layer 

Thickness, 

in. 

Modulus  of  Elasticity, 
psi 

Poisson's  Ratio 

Surface  Course 

3 

150,000 

0.25 

Base  Course 

6 

50,000 

0.30 

Subbase 

24 

25,000 

0.35 

Subgrade 

00 

5,000 

0.45 
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Item  4  of  the  flexible  pavement  test  section  was  analyzed  using  the  main 
program  of  the  AFPAV  code  and  the  theoretical  response  was  compared  with  the 
preliminary  experimental  test  data  obtained  at  WES  (ref.  4).  Table  I  presents 
the  layer  properties  that  were  used  for  the  various  layers  of  flexible  item  4 
in  the  analysis  reported  in  reference  4.  It  was  shown  therein  that  the  pave¬ 
ment  response  predicted  by  the  linear  elastic  analysis  did  not  vary  signifi¬ 
cantly  whether  item  4  was  considered  to  consist  of  seven  layers  as  shown  in  fig¬ 
ure  22a  or  to  consist  of  only  four  layers  as  idealized  in  table  I.  Therefore, 
item  3  of  the  flexible  test  section  was  idealized  as  a  four-layered  pavement 
system  and  analyzed.  The  layer  properties  of  item  3  are  shown  in  table  II. 

Since  the  value  of  108,000  psi  for  the  E-modulus  of  the  crushed-stone  base 


Table  I 

LAYER  PROPERTIES  OF  WES  FLEXIBLE  ITB4  4 


Layer 

Thickness , 
in. 

Modulus  of  Elasticity, 
psi 

Poisson's  Ratio 

Surface  Course 

3 

150,000 

0.25 

Base  Course 

6 

108,000 

0.25 

Subbase 

24 

22,500 

0.35 

Subgrade 

00 

4,400 

0.45 

Note:  Data  are  from  reference  4. 


Table  II 

ASSUMED  LAYER  PROPERTIES  OF  WES  FLEXIBLE  ITEM  3 


Layer 

Thickness , 
in. 

Modulus  of  Elasticity, 
psi 

Poisson's  Ratio 

Surface  Course 

3 

150,000 

0.25 

Base  Course 

6 

50,000 

0.30 

Subbase 

24 

25,000 

0.35 

Subgrade 

00 

5,000 

0.45 
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course  in  reference  4  was  too  high,  a  reduced  E-modulus  of  50,000  psi  was  as¬ 
sumed  for  this  layer  based  on  an  empirical  relationship  given  by  Seed,  et  al. 
(ref.  17).* 

Figure  23  shows  a  comparison  of  the  theoretical  stress  distribution  deter¬ 
mined  by  the  AFPAV  code  and  the  experimental  stress  distribution  measured  at 
WES  in  flexible  item  3  (ref.  15).  As  in  reference  4,  the  comparison  was  not 
satisfactory,  probably  because  of  the  inadequacy  of  the  linear  elastic  model 
currently  used  in  the  AITAV  code  and  the  choice  of  elastic  properties.  It  is 
hoped  that  by  using  the  strain -dependent  moduli  obtained  from  Hardin's  research 
(ref.  16)  for  the  base  and  subgrade  materials,  it  may  be  possible  to  achieve 
better  correlation  between  the  computed  and  the  measured  data;  the  AFPAV  code 
is  thus  being  modified  to  consider  strain -dependent  moduli. 

Figure  24  shows  a  comparison  of  the  theoretical  elastic  vertical  deflec¬ 
tions  and  the  experimental  data  reported  in  reference  15.  Although  the  agree¬ 
ment  between  computed  deflections  and  experimental  data  is  reasonably  good  at 
the  upper  layers  the  computed  deflection  curve  does  not  have  the  marked  change 
in  curvature  as  does  the  measured  deflection  curve.  Although  it  is  difficult 
at  this  stage  to  pinpoint  the  reasons  for  the  unsatisfactory  correlation  between 
the  measured  and  the  computed  data,  the  following  factors  may  be  responsible  to 
a  greater  or  lesser  degree  for  the  discrepancies: 

(1)  The  inability  to  define  the  most  appropriate  clastic  proper¬ 
ties  of  the  various  layers  of  the  pavement  system  (for  example, 
the  dependence  of  the  subgrade,  subbase,  and  base  course  ma¬ 
terials  on  the  actual  strain  level  is  not  considered  in  the 
linear  elastic  analysis  performed  in  this  investigation. — Xhc 
fact  that  the  granular  base  course  has  little  or  no  tensile 
strength  is  ignored  in  assuming  a  modulus  of  50,000  psi  in  all 
directions) ; 

(2)  The  subgrade  with  a  modulus  of  5,000  psi  is  assumed  to  extend 
to  a  depth  of  50  ft  from  the  surface  in  the  finite  clement 

f. 

Even  this  reduced  value  is  not  correct  since  it  should  be  zero  in  the  hori¬ 
zontal  direction  where  tensile  stress  is  developed  under  the  applied  loading. 
Since,  in  the  present  version  of  the  AFPAV  code,  different  moduli  values  can¬ 
not  be  input  in  different  directions,  predicted  response  will  be  in  error  to 
some  extent. 
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3 -in.  Asphalt 


Figure  23.  Comparison  of  Theoretical  Stress  Dis¬ 
tribution  and  WES  Field  Test  Data 


3- in.  Asphalt 

Maximum  Elastic  Vertical  Deflection,  in.  Cement  Con- 


Figure  24.  Comparison  of  Theoretical  Elastic  Vertical 
Deflections  and  WES  Field  Test  Data 
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model  while,  in  reality,  the  built-in  subgrade  with  a  CBR 
value  of  about  4  or  less  extends  to  a  depth  of  only  12  ft 
from  the  top  surface;  and 

(3)  The  actual  modulus  of  the  existing  ground  below  12  ft  is 
not  known,  but  it  is  probably  higher  than  that  of  the  CBR 
4  subgrade.  (If  so,  the  actual  deflection  will  be  less 
than  the  computed  values  at  greater  depths.) 

Suitable  modification  of  the  AFPAV  code  to  account  for  the  strain-dependent 
modulus  of  the  subgrade,  etc.,  and  no-tension  characteristic  of  the  granular 
material  is  being  made  now  to  improve  the  match  between  the  predicted  values 
and  the  measured  data. 

A  brief  clarification  is  required  to  explain  the  qualifying  term  maximum 
in  figures  23  and  24.  The  maximum  stress  or  deflection  means  the  maximum  quan¬ 
tity  at  a  given  depth  in  the  test  item  due  to  the  12-wheel  assembly  of  the  C-5A 
gear  used  in  the  WES  tests.  Each  wheel  carried  a  load  of  30,000  lb  with  a  con¬ 
tact  pressure  of  about  106  psi.  (The  stress  and  deflection  curves  shown  in 
figures  23  and  24,  respectively,  were  referred  to  as  limiting  curves  in  refer¬ 
ence  4.) 

2.  RIGID  PAVEMENT  TEST  SECTION 

Item  1  of  the  rigid  pavement  test  section  (fig.  22b)  was  analyzed  to  de¬ 
termine  its  response  to  static  loading  by  a  12-wheel  assembly  of  the  C-5A  gear. 

The  layer  properties  used  in  this  analysis  are  given  in  table  III. 

No  measured  vertical  stress  data  from  the  rigid  pavement  section  are  avail 

able  to  compare  with  the  computed  values;  however,  strain  data  were  received 


Table  III 

LAYER  PROPERTIES  OF  WES  RIGID  ITEM  1 


Layer 

Thickness, 

in. 

Modulus  of  Elasticity, 
psi 

Poisson's  Ratio 

Surface  Course 

10 

6,600,000 

0.20 

Subgrade  (CBR  4) 

590 

6,000 

0.45 

from  CERL  (ref.  18).  A  strain  of  23  yin. /in.  was  measured  at  the  center  of 
one  of  the  panels  of  item  1  when  it  was  statically  loaded  by  a  12 -wheel  assem¬ 
bly  of  the  C-5A  gear.  Each  wheel  had  a  load  of  15  kips.  One  of  the  two  inner 
wheels  of  the  interior  line  of  four  wheels  was  located  at  the  center  of  the 
slab  (ref.  4).  The  computed  strain  at  the  corresponding  location  was  30  yin./ 
in.  for  a  12-wheel  assembly  with  a  load  of  15  kips  in  each  wheel;  this  compares 
reasonably  well  with  the  measured  value. 

The  measured  recoverable  deformation  in  item  1  of  the  rigid  pavement  sec¬ 
tion  averaged  about  0.11  in.  at  the  start  of  traffic  tests  using  the  12-wheel 
assembly  of  the  C-5A  gear  at  WES  (ref.  15,  p.  27).  The  elastic  deflection  of 
item  1  under  the  same  loading  computed  by  the  AFPAV  code  was  also  0.11  in. 

This  shows  that,  when  the  pavement  system  follows  a  linear  elastic  law,  as  the 
layers  of  item  1  apparently  did  in  this  case,  good  correlation  between  theory 
and  actual  field  behavior  can  be  achieved  with  the  AFPAV  code.  This  statement 
is  qualified  by  saying  that  the  definition  of  the  appropriate  material  proper¬ 
ties  plays  a  vital  role  in  predicting  the  pavement  response. 


SUCTION  V 


PARAMETRIC  STUDY  OF  WES  TIiST  SHCTIONS 

To  demonstrate  the  fact  that  theoretical  predictions  of  pavement  response 
are  only  as  good  as  the  input  of  material  properties  and  fo  show  the  relative 
importance  of  the  various  parameters  of  pavement  systems,  a  parametric  study 
of  both  flexible  and  rigid  pavement  systems  was  conducted. 

1.  FLEXIBLE  PAVEMENTS 

The  five  items  of  the  full-scale  MWHGL  flexible  pavement  test  section  shown 
in  figure  22a  were  used  to  select  a  representative  list  of  pavement  parameters. 
Item  3  of  this  test  section,  whose  properties  arc  given  in  table  II,  was  chosen 
as  a  reference  system.  Table  IV  lists  the  parameter  variations  studied.  Study¬ 
ing  all  possible  combinations  of  the  variables  in  this  table  obviously  would  be 
time  consuming.  Therefore,  only  one  parameter  at  a  time  was  varied,  and  the 
change  in  the  pavement  response  due  to  that  variation  was  studied. 

In  this  parametric  study,  attention  was  focused  on  the  changes  noticed  in 
surface  deflections  and  in  significant  stresses  due  to  changes  in  the  proper¬ 
ties  of  the  pavement  components.  The  computed  pavement  response  due  to  the 
static  loading  of  a  12-wheel  assembly  of  the  C-5A  gear  with  30,000  lb  on  each 
wheel  was  used.  Tables  V,  VI,  and  VII  list  the  changes  in  deflections  and 
stresses  due  to  the  changes  in  the  pavement  parameters. 

r* 

a.  Effect  of  Young's  Modulus 

The  surface  clastic  deflections  were  not  changed  significantly  by  var¬ 
iations  in  the  E-modulus  of  the  surface  course,  base  course,  and  subbase  within 
the  range  of  values  shown  in  table  IV.  But,  variations  in  Young's  modulus  of 
the  subgrade  affected  the  surface  deflection  very  significantly.  For  example, 
a  67  percent  reduction  in  the  E-modulus  of  the  surface  course  fi.e.,  from 
1,500,000  to  500,000  psi)  resulted  in  an  increase  of  only  8  percent  in  surface 
deflection,  while  a  50  percent  reduction  in  the  subgrade  E-modulus  increased 
the  surface  deflection  by  about  60  percent.  This  is  because  in  layered  pave¬ 
ment  systems  the  softer  subgrade  contributes  more  to  the  surface  deflection  than 
the  stiffen  upper  layers  do  (ref.  12).  Figure  25  shows  the  effect  of  the  sub¬ 
grade  E-modulus  on  surface  deflections.  Even  though  increased  values  of  the  E- 
modulus  of  the  surface  course  reduced  the  vertical  stresses  in  the  base  course 
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Table  TV 


•LliXlBLU  PAVEMENT  PARAMETER  VARIATIONS  STUDIED 


Parameter  Varied 

Values  of  Parameter  Varied 

Symbol 

Unit 

E 

psi 

1.5  x  106 

500,000 

150,000 

V 

0.20 

0.25 

0.30 

h 

in. 

3.0 

4.5 

6.0 

E 

psi 

75,000 

50,000 

37,500 

V 

-- 

0.25 

0.30 

0.35 

h 

in. 

6.0 

9.0 

12.0 

E 

psi 

37,500 

25,000 

12,500 

\) 

-- 

0.30 

0.35 

0.40 

h 

in. 

15.0 

24.0 

33.0 

E 

psi 

5,000 

2,500 

1,250 

V 

-  - 

0.40 

0.45 

0.495 

E 

psi 

5,000 

50,000 

500,000 

Layer 


Surface  Course 
Surface  Course 
Surface  Course 
Base  Course 
Base  Course 
Base  Course 
Subbase 
Subbasc 
Subbasc 
Subgrade 
Subgrade 
Lxisting  Ground 
Profile 


Note:  h  =  layer  thickness. 


Table  V 

EFFECT  OF  PARAMETER  VARIATIONS  IN  SURFACE  COURSE 
ON  FLEXIBLE  PAVEMENT  RESPONSE 


Parameters 

Pavement  Responses 

Remarks 

■ 

V 

Maximun  Elastic 
Surface  Deflection 
in. 

Maximum  Hori¬ 
zontal  Stress, 
psi 

Compression 

Tension 

1,500,000 

500,000 

150,000 

150,000 

150,000 

150,000 

150,000 

150,000 

150,000 

0.25 

0.25 

0.25 

0.20 

0.25 

0.30 

0.25 

0.25 

0.25 

0.157 

0.170 

0.183 

0.183 

0.183 

0.182 

0.183 

0.174 

0.168 

435 

279 

166 

155 

166 

179 

166 

150 

136 

191 

23 

0 

0 

0 

0 

0 

0 

0 

E -modulus  only 
varied 

Poisson's  ratio 
only  varied 

Thickness  only 
varied 

Note:  The  properties  of  the  other  layers  are  given  in  table  II. 


Table  VI 


hi' FI iCl  OF  l'ARAMimiR  VARIATIONS  IN  BASF  COURSF 
AND  SUBBASF  ON  FLEXIBLE  PAVB1ENT  RESPONSE 


Parameters 

Pavement  Responses 

E, 

V 

b, 

Maximum  Elastic 

Maximum  Vertical 

Remarks 

Surface  Deflection, 

Stress 

in  Subgrade, 

psi 

in. 

in. 

psi 

Base  Course 

37,500 

50,000 

0.30 

0.30 

6 

6 

0.187 

0.183 

8.7 

8.5 

E-modulus  only 
varied 

75,000 

0.30 

6 

0.178 

8.4 

50,000 

50,000 

50,000 

0.25 

0.30 

0.35 

6 

6 

6 

0.183 

0.183 

0.182 

8.5 

8.5 

8.5 

Poisson's  ratio 
only  varied 

50,000 

50,000 

0.30 

0.30 

6 

9 

0.183 

0.174 

8.5 

7.7 

Thickness  only 
varied 

50,000 

0.30 

12 

0.167 

7.1 

Subbase 

12,500 

25,000 

0.35 

0.35 

24 

24 

0.207 

0.183 

9.5 

8.5 

E -modulus  only 
varied 

37,500 

0.35 

24 

0.171 

8.0 

25,000 

25,000 

25,000 

0.30 

0.35 

0.40 

24 

24 

24 

0.182 

0.183 

0.183 

8.5 

8.5 

8.6 

Poisson's  ratio 
only  varied 

25,000 

25,000 

25,000 

0.35 

0.35 

0.35 

15 

24 

33 

0.197 

0.183 

0.168 

- 

11.7 

8.5 

6.7 

Thickness  only 
varied 

Note:  The  properties  of  the  other  layers  are  given  in  table  II. 


Table  VII 


EFFECT  OF  PARAMETER  VARIATIONS  IN  SUBGRADE 
ON  FLEXIBLE  PAVEMENT  RESPONSE 


Parameters 

Maximum  Elastic 

li, 

V 

Surface  Deflection, 

Remarks 

psi 

in. 

1 

BE  ; 

BJfggjM 

0.450 

0.450 

0.450 

0.183 

0.299 

0.493 

E-modulus  only 
varied 

1 

Mrt? ' , 

* 

0.215 

0.183 

0.141 

Poisson's  ratio 
only  varied 

Note:  The  properties  of  the  other  layers  are  given  in  table  II. 


and  subbase,  they  did  not  alter  the  subgrade  stresses  significantly  (fig.  26). 
This  also  explains  why  changes  in  the  surface  course  E-modulus  within  the 
ranges  shown  in  table  IV  do  not  affect  the  surface  deflection  significantly. 
When  the  E-modulus  was  increased  tenfold  (from  150,000  psi  to  1,500,000  psi), 
the  decrease  in  surface  deflection  was  only  about  14  percent  (table  V). 

The  horizontal  stress  in  the  surface  course  was  affected  by  its 
Young's  modulus  as  shown  in  figure  27.  With  an  E-modulus  of  500,000  psi,  a 
horizontal  tensile  stress  of  about  100  psi  was  developed,  while  with  an  E- 
modulus  of  1,500,000  psi,  a  larger  horizontal  tensile  stress  (greater  than  200 
psi),  was  developed  because  of  increasing  beam  action.  There  was  no  horizon¬ 
tal  tensile  stress  in  the  surface  course  when  the  E-modulus  of  the  layer  was 
150,000  psi.  But  this  was  true  only  when  the  E-modulus  of  the  base  course  was 
50,000  psi  and  the  E-modulus  of  the  subbase  was  25,000  psi.  An  E-modulus  of 
150,000  psi  for  the  surface  course,  with  a  base  course  F.-modulus  of  25,000  psi 
and  a  subbase  E-modulus  of  15,000  psi,  produced  a  horizontal  tensile  stress  of 
about  30  psi  in  the  surface  course  (ref.  4).  This  shows  that  the  modular  ra¬ 
tios  (i.e.,  ratios  of  Young's  moduli)  of  the  various  pavement  components  are 
important  parameters  which  influence  the  pavement  response.  The  modular  ratio 
of  the  surface  course  to  the  base  course  affects  the  shear  stresses  in  these 
layers  significantly  (fig.  28). 
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Figure  26.  Effect  of  Surface  Course  E-Modulus  on  Vertical  Stress 
in  Flexible  Pavement 


Compression  Horizontal  Stress,  psi  Tension 


Effect  of  Surface  Course  E-Modulus  on  Horizontal  Stress 
in  Flexible  Pavement 


ea 


gure  28. 


Even  though  the  E-modulus  values  of  the  base  course  and  subbase  may 
not  be  very  significant  for  surface  elastic  deflections,  they  are  very  impor¬ 
tant  in  determining  the  state  of  stress  in  the  upper  layers  of  the  pavement 
system.  Therefore,  it  is  necessary  to  determine  the  E-moduli  of  the  various 
layers  as  accurately  as  possible  since  structural  failure  of  the  pavement  sys¬ 
tem  can  be  caused  by  critical  load  stresses  exceeding  the  inherent  strength  of 
the  pavement  materials. 

b.  Effect  of  Poisson's  Ratio 

Surface  deflections  were  not  significantly  affected  by  Poisson's  ra¬ 
tio  of  the  surface  course,  base  course,  or  subbase  for  the  range  of  thicknesses 
studied;  however,  Poisson's  ratio  of  the  subgrade  influenced  the  deflections 
significantly  (fig.  29).  A  Poisson's  ratio  of  0.40  for  the  subgrade  gave  a 
maximum  surface  deflection  of  about  0.22  in.,  while  a  Poisson's  ratio  of  0.495 
gave  a  maximum  deflection  of  only  0.14  in.  (These  maximum  values  are  shown  in 
table  VII). 

c.  Effect  of  Layer  Thickness 

Increasing  the  thicknesses  of  the  upper  pavement  layers  (i.e.,  the 
surface  course,  base  course,  and  subbase)  reduced  the  vertical  stresses  in  the 
subgradc  and,  consequently,  decreased  the  surface  deflections.  For  example, 
when  the  subbase  thickness  (originally  24  in.)  was  altered  to  15  and  33  in., 
the  vertical  stress  at  the  top  of  the  subgrade  became  11.7  and  6.7  psi,  respec¬ 
tively  (fig.  30  and  table  VI).  This  results  in  a  reduction  of  about  15  percent 
in  surface  deflection.  The  maximum  elastic  surface  deflection  was  about  0.20 
in.  with  a  subbase  thickness  of  33  in.,  while  it  was  0.17  in.  with  a  subbasc 
thickness  of  15  in.  (table  VI). 

d.  Effect  of  Young's  Modulus  of  Existing  Ground 

In  the  finite  element  idealization  of  the  pavement  system,  the  total 
depth  of  the  pavement  structure  was  assumed  to  be  50  ft  and  the  subgradc  was 
considered  to  have  a  constant  E-modulus.  The  total  built-up  depth  of  the  WES 
flexible  test  section  (fig.  22a)  was  only  12  ft.  The  E -modulus  of  the  ground 
below  the  built-up  subgrade  was  not  known.  However,  based  on  the  reported  CBR 
value  of  2  to  4,  the  E-modulus  of  the  ground  was  assumed  to  be  the  same  as  that 
of  the  subgrade.  To  illustrate  the  importance  of  knowing  the  E-modulus  of  the 
existing  ground  below  a  depth  of  12  ft,  a  parametric  study  assuming  three 
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^  _ ^  (c)  v  =  0.495  "* 

Note:  +  indicates  the  wheel  positions  of  the 
12-wheel  assembly  of  the  C-5A  gear. 

Effect  of  Subgrade  Poisson's  Ratio  on  Surface 
Deflection  in  Flexible  Pavement 


Figure  29 


trriVBl 

LuiwJ 

MiIMtvi 

Vertical  Stress,  psi 
10  20  30  40  50 


different  moduli  values  (5,000,  50,000,  and  500,000  psi)  was  conducted.  The 
maximum  surface  deflection  of  0.185  in.  with  an  H-modulus  of  5,000  psi  was  re¬ 
duced  by  one-third  to  0.121  in.  when  the  existing  ground  was  assumed  to  have 
an  li -modulus  of  50,000  psi.  An  E-modulus  of  500,000  psi  below  12  ft  further 
reduced  the  maximum  surface  deflection  to  0.108  in.,  a  reduction  of  about  40 
percent  (fig.  31).  Therefore,  it  appears  that  it  is  necessary  to  determine  the 
properties  of  the  pavement  system  to  an  appreciable  depth. 

2.  RIGID  PAVEMENTS 

Table  VIII  summarizes  the  various  sets  of  layer  properties  used  in  the 
parametric  study  of  rigid  pavement  systems.  The  computed  pavement  response 
due  to  the  static  loading  of  a  12-wheel  assembly  of  the  C-5A  gear  with  30,000 
lb  on  each  wheel  was  used.  The  effect  of  the  parameter  variations  on  two  im¬ 
portant  responses,  viz.,  the  maximum  surface  deflections  and  the  maximum  ten¬ 
sile  stresses  in  the  surface  course  are  listed  in  tables  IX,  X,  and  XI. 

a.  Effect  of  Young's  Modulus  of  Surface  Course 

Figure  32  shows  the  effect  of  varying  the  E-modulus  of  the  surface 
course  on  the  deflection  basin  of  the  pavement  structure  caused  by  the  static 
loading  of  the  12-wheel  assembly  of  the  C-5A  gear.  Decreasing  the  E-modulus 
of  the  surface  course  from  6,600,000  to  3,000,000  psi  produced  about  an  18- 
percent  increase  in  the  maximum  surface  deflection. 

Figure  33  shows  how  varying  the  E-modulus  of  the  surface  course  af¬ 
fects  the  distribution  of  vertical  normal  stress,  horizontal  normal  stress, 
and  horizontal  shear  stress.  In  this  analysis,  the  E-modulus  of  the  subgrade 
was  kept  constant  at  6,000  psi,  and  the  thickness  of  the  surface  course  was 
held  constant  at  10  in.  Increasing  the  surface  course  E-modulus  from  3,000,000 
to  6,600,000  psi  increased  the  maximum  horizontal  tensile  stress  in  the  sur¬ 
face  course  by  about  25  percent  (table  IX).  This  increasing  plate  action  of 
the  surface  course  decreases  the  maximum  vertical  stress  in  the  subgrade.  No 
change  in  shear  stress  was  noticed  for  the  range  of  surface  course  E-moduli 
studied. 

b.  Effect  of  Layer  Thickness 

Increasing  the  thickness  of  the  surface  course  from  10  to  12  and  14 
in.  decreased  the  maximum  deflection  by  about  17  percent  (table  IX  and  fig.  34). 
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Table  VI 11 


RIGID  PAVEMENT  PARAMETER  VARIATIONS  STUDIED 


Surface  Course  Parameters 

Subgrade  Parameters 

Thick- 

Modulus  of 

Poisson's 

Layer 

Thick- 

Modulus  of 

Poisson's 

ness  (h) , 

Elasticity  (E) , 

Ratio  (v) 

ness  (h) , 

Elasticity  (E) , 

Ratio  (v) 

in. 

psi 

in. 

psi 

10.0 

6.6  x  106 

0.20 

1 

600 

6,000 

0.45 

10.0 

5  x  106 

0.20 

1 

600 

6,000 

0.45 

10.0 

3  x  106 

0.20 

1 

600 

6,000 

0.45 

10.0 

6.6  x  106 

0.20 

1 

600 

6,000 

0.45 

12.0 

6.6  x  106 

0.20 

1 

600 

6,000 

0.45 

14.0 

6.6  x  106 

0.20 

1 

600 

6,000 

0.45 

10.0 

6.6  x  106 

0.20 

1 

600 

6,000 

0.45 

10.0 

6.6  x  106 

0.20 

1 

600 

3,000 

0.45 

10.0 

6.6  x  106 

0.20 

1 

600 

1,000 

0.45 

10.0 

6.6  x  106 

0.20 

1 

2 

12 

588 

24,000 

6,000 

0.35 

0.45 

10.0 

6.6  x  106 

0.20 

1 

2 

18 

582 

36,000 

6,000 

0.35 

0.45 

10.0 

6.6  x  106 

0.20 

1 

2 

24 

576 

60,000 

6,000 

0.35 

0.45 

Table  IX 

DEFECT  OF  PARAMETER  VARIATIONS  IN  SURFACE  COURSE 
ON  RIGID  PAVEMENT'  RESPONSE 


Parameters 

Pavement  Responses 

E, 

h, 

Maximum  Elastic 

Maximum  Horizontal 

Remarks 

Surface  Deflection, 

Tensile  Stress, 

psi 

in. 

in. 

psi 

6.6  x  106 
5.0  x  106 

10 

10 

0.098 

0.105 

418 

388 

L-modulus  only 
varied 

3.0  x  106 

10 

0.116 

335 

6.6  x  106 
6.6  x  106 
6.6  x  10 6 

10 

i: 

14 

0.098 

0.088 

0.081 

418 

338 

278 

Thickness  only 
varied 

Note:  The  surface  course  rested  directly  on  the  subgrade;  the  E-modulus  and 
Poisson's  ratio  of  the  subgrade  were  6000  psi  and  0.45,  respectively. 
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Table  X 


EFFECT  OF  PARAMETER  VARIATIONS  IN  SUBGRADE 
ON  RIGID  PAV04ENT  RESPONSE 


Parameters 

Pavement  Responses 

E, 

psi 

V 

Maximun  Elastic 
Surface  Deflection, 
in. 

Maximum  Horizontal 
Tensile  Stress, 
psi 

6,000 

0.45 

0.098 

418 

3,000 

0.45 

0.166 

488 

1,000 

0.45 

0.387 

S79 

Note: 

The  surface  course 
on  the  subgrade. 

with  E  *  6.6  x  106  psi  and 

v  =  0.20  rested  directly 

Table  XI 

EFFECT  OF  BASE  COURSE  ON  RIGID  PAVB4ENT  RESPONSE 


Base  Course  Parameters 

Pavement  Responses 

h. 

E, 

Maximum  Elastic 

Maximum  Horizontal 

V 

Surface  Deflection, 

Tensile  Stress, 

in. 

psi 

in. 

psi 

12 

24,000 

0.35 

0.093 

364 

18 

36,000 

0.35 

0.090 

341 

24 

60,000 

0.35 

0.084 

298 

Note:  The  surface  course  parameters  were  E  -  6.6  x  10 6  psi,  v  *  0.20,  and 

h  =  10  in.;  the  subgrade  parameters  were  E  ■  6  x  103  psi,  and  v  *  0.45. 


In  these  two  cases  the  properties  of  the  subgrade  as  well  as  the  elastic  con¬ 
stants  of  the  surface  course  remained  unchanged. 

r 

Figure  35  shows  the  stress  distribution  in  surface  courses  of  varying 
thicknesses  (10,  12,  and  14  in.).  The  E-modulus  of  the  surface  course  and  that 
of  the  subgrade  were  held  constant.  The  magnitude  of  maximum  horizontal  ten¬ 
sile  stress  in  the  14-in. -thick  surface  course  was  only  about  two-thirds  of 
that  in  the  10-in. -thick  surface  course.  Similarly,  the  vertical  normal  stress 
at  the  14 -in. -thick  surface  course/subgrade  interface  was  also  reduced  by  about 
25  percent  (fig.  35)  from  that  of  the  10-in. -thick  surface  course,  thus  result¬ 
ing  in  decreased  deflection  (table  IX) . 
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Horizontal  Normal  Stress,  psi  Vertical  Normal  Stress,  psi  Horizontal  Shear  Stress,  psi 
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Effect  of  Surface  Course  E-Modulus  on  Stress  Distribution  in  Rigid  Pavement 


Note:  +  indicates  the  wheel  positions  of  the  12-wheel 
assembly  of  the  C-5A  gear. 

Figure  34.  Effect  of  Surface  Course  Thickness  on 
Surface  Deflection  in  Rigid  Pavement 


Vertical  Normal  Stress 


c.  Effect  of  Young's  Modulus  of  Subgrade 


f  ig-ire  36  shows  the  effect  of  changing  the  E-modulus  of  the  subgradc 
while  keeping  the  surface  course  properties  constant.  Reducing  the  subgrade 
stiffness  from  6,000  to  3,000  psi  increased  the  maximum  deflection  by  about  70 
percent,  while  reduction  from  6,000  to  1,000  psi  increased  the  maximum  deflec¬ 
tion  by  about  400  percent  (table  X) . 

Figure  37  shows  the  changes  in  stress  distribution  caused  by  varia¬ 
tion  of  subgrade  E-modulus  from  1,000  to  6,000  psi.  The  weaker  the  sul  Tide 
the  greater  the  horizontal  normal  stress  in  the  surface  course.  The  maximum 
horizontal  tensile  stress  in  the  surface  course  was  increased  by  about  40  per¬ 
cent  when  the  subgradc  E-modulus  was  reduced  from  6,000  to  1,000  psi  (tabic  X). 
The  vertical  normal  stress  at  the  surface  course/subgrade  interface  was  reduced 
by  about  25  percent  for  the  same  reduction  in  subgrade  E-modulus.  However,  in 
spite  of  this  reduction  in  subgrade  stress,  the  surface  deflection  was  increased 
by  about  400  percent  because  of  the  reduction  in  subgrade  E-modulus  from  6,000 
to  1,000  psi. 

d.  Effect  of  Base  Course 

Figure  38  and  table  XI  show  the  manner  in  which  the  deflection  basin 
of  the  surface  course  is  affected  when  there  is  a  base  course  below  the  surface 
course.  While  the  deflection  basin  of  the  surface  course  resting  on  a  12-in.- 
thick  base  course  with  an  E-modulus  of  24,000  psi  varied  only  slightly  from 
that  of  a  similar  surface  course  resting  directly  on  the  subgrade,  a  24-in. - 
thick  base  course  with  an  E-modulus  of  60,000  psi  reduced  the  maximum  deflec¬ 
tion  by  about  10  percent  and  made  the  deflection  basin  somewhat  flatter.  In 
both  cases,  the  E-modulus  of  the  surface  course  was  6,600,000  psi  and  that  of 
the  subgrade  was  6,000  psi. 
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Figure  37.  Effect  of  Subgrade  E-Modulus  on  Stress  Distribution  in  Rigid  Pavement 


SECTION  VI 


PARAMETRIC  STUDY  FOR  AIR  FORCE  CIVIL  ENG  INFER  INC  CENTER 

The  IVES  parametric  study  of  both  flexible  pavement  and  rigid  pavement  sys- 

* 

terns  (section  V)  was  conducted  from  a  theoretical  point  of  view,  assuming  some 
representative  values  of  the  pavement  layer  properties.  In  addition  to  this 
study,  which  was  reqi  ired  by  this  research  effort,  another  parametric  study  of 
a  rigid  pavement  system  was  conducted  at  the  request  of  the  Civil  Engineering 
Center  (CEC)  at  Wright -Patterson  AFB.  This  afforded  an  opportunity  to  demon¬ 
strate  the  practical  use  of  the  AFPAV  code. 

Figure  39  shows  the  pavement  system  parameters  used  in  the  study  for  the 
CEC.  The  objective  of  the  study  was  the  determination  of  the  maximum  horizon¬ 
tal  tensile  stress  in  the  surface  course  when  loaded  by  (1)  the  C-5A,  (2)  the 
Boeing  747,  and  f3)  the  C-141A  aircraft.  Figure  40  shows  the  configuration  of 
the  main  landing  gears  of  the  C-141A  and  Boeing  747  aircraft;  table  XU  com¬ 
pares  the  assumed  wheel  data  for  the  three  aircraft. 

Figure  41  shows  the  effect  of  the  surface  course  thickness  on  the  maximum 
horizontal  tensile  stress  in  the  rigid  pavement  system.  The  tensile  stresses 
caused  by  the  C-5A  aircraft  were  lower  than  those  caused  by  the  other  two  air¬ 
craft  because  the  arrangement  of  the  wheels  of  the  landing  gears  (fig.  5) 

and  the  low  tire  pressure  (table  XII)  of  the  C-5A  aircraft.  Although  for  thick¬ 
nesses  less  thar  12  in.,  the  C-141A  produced  larger  tensile  stresses  than  did 
the  Boeing  747  aircraft,  this  situation  reversed  when  the  thickness  of  the  sur¬ 
face  course  was  increased.  In  this  analysis,  as  in  most  of  the  analyses  in¬ 
cluded  in  this  report,  it  was  assumed  that  the  loading  was  in  the  interior  por¬ 
tion  of  a  wide  pavement  system.  Furthermore,  the  presence  of  joints  in  the 
surface  course  was  also  ignored. 


The  results  of  this  study  have  been  published  (refs.  19,20). 
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Figure  39.  Rigid  Pavement  System  Analyzed  for  CEC 
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Table  XII 


DATA  FOR  MAIN  LANDING  GEARS  IN  CEC  STUDY 


Aircraft 

Number  of 
Wheels 

Contact  Area,* 
in.2 

Contact  Pressure,* 
psi 

Wheel  Load,* 
kips 

C-5A 

24 

285 

106 

Boeing  747 

16 

203 

204 

C-141A 

8 

208 

180 

37.4 

* 

Assumed  data. 
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Note:  d  is  the  depth  at  which  stresses 
shown  in  this  figure  existed. 

Figure  41.  Effect  of  Surface  Course  Thickness  on  Maximum  Horizontal 
Tensile  Stress  in  Surface  Course  of  Rigid  Pavement 
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SECTION  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 


1.  CONCLUSIONS 

Subject  to  the  obvious  limitations  inherent  in  a  linear  elastic  model  of 
layered  pavement  systems,  the  AFPAV  code  appears  to  be  a  satisfactory  analyti¬ 
cal  tool  for  predicting  the  pavement  response  due  to  the  MWI1GL  of  modem  air¬ 
craft  such  as  the  C-5A.  This  is  supported  by  the  results  of  the  analysis  of 
the  full-scale  pavement  sections  (particularly  the  rigid  pavement  section) 
tested  at  WES. 

Even  though  accurate  predictions  of  pavement  response  may  be  difficult  to 
make  at  present,  either  because  of  the  limitations  of  the  code  or  because  of 
the  engineer's  inability  to  determine  the  proper  input  parameters  pertaining 
to  the  different  pavement  layers ,  the  AFPAV  code  is  a  very  useful  tool  for 
parametric  studies  of  pavement  systems  loaded  by  multiple  wheels.  The  results 
of  the  parametric  studies  of  both  rigid  and  flexible  pavements  indicate  the 
need  for  accurate  determination  of  the  following  pavement  system  parameters  for 
a  linear  elastic  model: 

(1)  Moduli  of  elasticity  of  all  pavement  layers 

(2)  Poisson's  ratio  of  the  subgrade 

(.3)  Thicknesses  of  the  upper  layers 

(4)  Elastic  constants  of  the  existing  ground  below 
the  built-up  sub grade 

2.  RECOMMENDATIONS 

Experience  in  working  with  the  AFPAV  code  indicates  the  need  for  the  fol¬ 
lowing  improvements  in  the  code: 

(1)  The  code  should  be  modified  to  consider  the  strain -dependent  moduli 
values  of  the  pavement  layers. 

f 2 )  The  phenomenon  of  no-tens'  m  in  certain  layers  should  be  considered. 
(In  such  an  analysis,  the  effects  of  the  geostatic  stresses  (ref.  21)  should 
be  included.) 
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(3)  The  code  should  be  modified  to  analyze  the  sliding -plane  problem, 
both  along  vertical  joints  and  along  interfaces,  when  the  pavement  is  loaded 
by  multiple  wheels. 

(4)  Viscoelastic  behavior  of  asphaltic  concrete  layers  as  well  as  temper¬ 
ature  effects  should  be  analyzed. 

(5)  To  predict  permanent  deformations,  the  unloading  phenomenon  should  be 
analyzed.  (So  far  attention  has  been  paid  only  to  the  loading  portion  of  the 
stress/strain  relationships.  This  gives  only  the  instantaneous,  recoverable  de¬ 
formations.) 

(6)  Dynamic  analyses  of  pavement  response  using  the  AFPAV  code  should  be 

undertaken  to  see  if  there  are  any  conditions  which  might  be  more  critical  than 

* 

the  static  analyses.  (This  may  eventually  lead  to  the  refinement  of  analyti¬ 
cal  studies  of  runway  roughness  which  at  present  do  not  consider  the  layered 
pavement  systems.) 

Modification  of  the  AFPAV  code  to  perform  the  analyses  in  (1),  (2),  and 
(3)  is  currently  in  progress.  When  the  modified  code  becomes  available,  it 
should  be  possiole  to  predict  the  three  basic  types  of  pavement  response  (i.e., 
deflection,  strain,  and  stress)  more  accurately.  Thus,  the  AFPAV  code  will 
be  of  great  help  in  formulating  workable  pavement  distress  criteria  for  pave¬ 
ment  evaluation. 


■k 

Lockheed-Califomia  Company  has  performed  some  research  for  the  Federal  Avia¬ 
tion  Administration  to  predict  the  dynamic  wheel  load  effects  on  airport  pave¬ 
ments  using  the  plate  on  elastic  solid  theory  (ref.  22). 
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